Substrate concentration dependence of the diffusion-controlled steady-state rate 
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The Smoluchowski approach to diffusion-controlled reactions is generalized to interacting sub- 
strate particles by including the osmotic pressure and hydrodynamic interactions of the nonideal 
particles in the Smoluchoswki equation within a local-density approximation. By solving the strictly 
linearized equation for the time-independent case with absorbing boundary conditions, we present 
an analytic expression for the diffusion-limited steady-state rate constant for small substrate con- 
centrations in terms of an effective second virial coefficient • Comparisons to Brownian dynamics 
simulations excluding HI show excellent agreement up to bulk number densities of BJpo ^ 0.4 for 
hard sphere and repulsive Yukawa-like interactions between the substrates. Our study provides an 
alternative way to determine the second virial coefflcient of interacting macromolecules experimen- 
tally by measuring their steady-state rate constant in diffusion-controlled reactions at low densities. 



I. INTRODUCTION 

Simple irreversible reactions of the type A+B^ A 
(with A the sink and B the substrates) are commonly 
found in many (bio)chemical processes, such as fluores- 
cence quenching, enzyme catalysis, polymerization, or 
colloid and protein aggregation, just to mention a few 
examplesfl The key parameter for these processes is the 
reaction rate constant, a measure for the number of re- 
actions per unit time, first addressed in the pioneering 
and now classical works of SmoluchowskiSi and Debye^ 
decades ago. Since then, various improvements and re- 
finements have been made in predicting the rate con- 
stants for diverse reactionS) ^1'^1^^^'^'^'^'^"d^1^^ in partic- 
ular, including solvent-mediated hydrodynamic interac- 
tions (HI) between sink (A) and substrate (B)(i or ex- 
amining effects of a non-zero concentration of the sink 
particlesi ^i^^i^"^ In most of the previous studies, interac- 
tions between the substrate particles were ignored, which 
is only justified in the case of very weakly interacting sub- 
strate particles or at infinite dilution. In recent attempts 
the infiuence of the excluded volume of the substrate par- 
ticles was examinedi^si^ and predicted to yield an in- 
creased reaction rate with increasing excluded volume or 
substrate concentration. Along these lines, Senapati and 
McCammonii also gave evidence for a strong influence of 
substrate interactions (without HI) on the reaction rate 
by means of Brownian Dynamics (BD) computer simula- 
tions. We tie up to these studies in this work and aim at 
a systematic examination of the substrate concentration- 
dependence of the rate constant, while the sinks remain 
at infinite dilution. 

For many of the reactions mentioned above the rate- 
limiting factor is the diffusional encounter of the reac- 
tants, particularly when the subsequent transformation 
does not involve a large activation barrier. In our work 
we will focus on the steady-state (long-time) rate of these 
diffusion-limited or diffusion-controlled reactions. The 
usual theoretical approach is based on the Smoluchowski 
equation (SE) which becomes time-independent in the 



steady-state. For weakly interacting particles or for small 
concentrations, the inhomogeneous density profile of the 
substrates around the sink varies smoothly over distances 
larger compared to the typical interaction range and well 
established generalizations of the SE are availableji^ of- 
ten used for instance for the problem of colloidal sedimen- 
tation. The generalized SE employs the osmotic pressure 
and density-dependent mobility of the interacting parti- 
cles within a local-density approximation (LDA)iS which 
assumes local homogeneity of the density and is justi- 
fied whenever the density varies slowly in space. Using a 
simple model with a spherical, isotropically reactive sink 
particle, we will show that a strictly linearized version of 
the generalized SE allows us to write down an analytic 
solution for the first order correction of the rate constant 
in linear order of substrate concentration. As a result, 
the correction coefficient is basically given by the second 
virial coefficient i?2, the first order correction coefficient 
in the expansion of the virial equation of state, but must 
be corrected for HI and is also infiuenced by the inter- 
action between sink and substrate. Fortunately, most 
chemical reactions occur at small densities of the reac- 
tants so that our result should be valid for a wide range 
of processes and systems. Despite the simplicity of our 
model the results are general and should be applicable to 
more realistic and complicated systems. 

In order to examine the range of validity of our theo- 
retical result we perform Brownian dynamics (BD) simu- 
lations for different systems in which the interaction be- 
tween the B particles is varied from hard sphere to at- 
tractive and repulsive Yukawa-like interactions, the lat- 
ter typically found in ionic solutions. The BD simula- 
tions do not include HI, which can be accurately treated 
only by more sophisticated and computationally more ex- 
pensive means, such as Lattice-Boltzmann methods-- or 
others iSiiSSiS^ However, in comparison to the theory ex- 
cluding HI we find excellent agreement up to substrate 
densities i?2Po — 0.4, showing the reliability of our theo- 
retical concept, and proving its applicability to the case 
when HI can be neglected, i.e. for long-ranged inter- 
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actions. In principle, our study opens up an alterna- 
tive way to determine the second virial coefficient of pro- 
teins or other particles by measuring the reaction rate in 
diffusion-limited reactions at small substrate densities. 

The paper is organized as follows: In section II we 
present the basic equations of motion and approxima- 
tions of our theory and arrive at a first order linear dif- 
ferential equation which can be solved analytically. The 
solution for the steady-state rate constant and density 
profile are presented and discussed in section III. A sys- 
tematic comparison to BD simulations in order to exam- 
ine the range of validity of our theory follows in section 
IV. Sec. V concludes our work with a few final remarks. 



II. GENERALIZED SMOLUCHOWSKI THEORY 

A. Model and basic equations 

Let us consider a spherical, isotropically reactive sink 
particle A with diameter ctaa which reacts with a B par- 
ticle (substrate) of size ctbh ='■ c when they touch at a 
center-to-center distance r = ctab = (caa + o')/2. The 
sinks are at infinite dilution while the B particles have a 
bulk number density po- Furthermore a potential Vab('') 
is acting between the sink and the substrate particles. 
Both A and B are dispersed in a solvent, which is taken 
into account by a position-dependent diffusion constant 
-Do(r) of the substrates, reflecting the diffusion of an iso- 
lated B particle relative to the A particle. We assume 
the sink to be completely at rest, which is justified for 
large sinks ctaa ^ cr. For close distances to A, Do{r) is 
determined by the HI between A and B, while for large 
distances when HI are negligible, we obtain by the 
Stokes-Einstein relation 



Doi\r\ -^^)=Do = (3/37r77CT)-i, 



(1) 



where = kgT is the thermal energy, and rj the sol- 
vent viscosity. 

The time-dependent Smoluchowski equation (SE) for 
noninteracting particles moving in a potential VAB(r) 
with a position-dependent diffusion constant Do{r) is 
given by 



dp{r,t)/dt = -V-J(r,t) 
= V-i^o(r) 



(2) 



t) — + Vp(r, t) 



For weak density inhomogeneities, which is always true 
for small densities or weak interactions, one can introduce 
substrate interactions by using the generalized Stokes- 
Einstein relation^'* 



D{p) = M{p) 



dH(p) 
dp ■ 



(3) 



where D{p) is the collective diffusion coefficient, M{p) is 
the density dependent mobility of the substrate particles, 



and n(p) the osmotic pressure of the interacting parti- 
cles. This is a generalized Stokes-Einstein equation in 
the sense that it generalizes (1) to the case of interacting 
particles in a homogeneous solution. The mobility M{p) 
is a reciprocal friction and is defined as the proportional- 
ity constant between the drift velocity and total force on 
a Brownian particle [velocity = M{p) x total force) 
in a steady-state situation. Within the LDA the gen- 
eralized Stokes-Einstein relation is applied to the local, 
position-dependent density p(r) of the substrate parti- 
cles, which we assume to vary smootly over distances 
larger compared to the range of the interaction poten- 
tial. The equation of motion can now be written as^ 



(9p(r,t)/(9i = - V- J(r,t) 



(4) 



= V-Af(p(r,<),r) 



P(r,t)- 



dr 



[Vp(r,t)] 



dH(p(r,t) 
dp{v) 



and is a generalization of Eq. ^ to weakly interacting 
substrate particles. The last term on the rhs of Eq. Q 
accounts for the force on a Brownian particle due to an 
unbalanced osmotic pressure caused by a concentration 
gradient of interacting particles. In addition, the density- 
dependent mobility corrects for HI between B particles. 
Note that for ideal particles dH/dp — ksT, and with no 
HI between B particles M(p(r), r) = (3Dq{y) and we find 
(2) again. To account for reactions at the sink we have 
to solve Eq. with the boundary condition 



(5) 



where ki is the intrinsic rate constant. In the limit of 
fully absorbing boundary conditions, ki — > oo, while 
p(o'ab) — > 0, and the reaction rate is only limited by 
the diffusion of the approaching substrates. We are in- 
terested in the steady-state solution for the spherically 
symmetric problem, with which the equation of motion 
reduces to 







d_ 

dr 



r^M{p{r),r) 



or or dp(rj 



with the boundary conditions for diffusion-controlled re- 
actions 



and 



47rCTABjVAB) = kpo, 



p(ctab) = and p(oo) = po, 



(7) 



(8) 



where j(ctab) is the absolute flux perpendicular through 
the sink surface at radial distance ctab from the sink cen- 
ter, and the quantity k is the diffusion-controlled steady- 
state rate constant which we are interested in. 

The osmotic pressure can be expressed in terms of the 
equation of state, or in general, the virial equation, which 
is an expansion of the pressure in terms of the particle 
density p. For small densities the series is usually trun- 
cated after the second order and the osmotic pressure 
reads 



(6) 



f3U{p) = P + B2P^ 



(9) 
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Conveniently, the second virial coefficient B2 can be cal- 
culated from the pair interaction VBB(r) =: ^(r) between 
the substrates viaASi 



Bo 



dr[l - exp(-/3y(r))]. 



(10) 



The local mobility of the substrates will be in general 
reduced by a nonzero concentration of other surrounding 
particles due to HI, and for small densities the correction 
is linear in p, so that we can write 



M{p{r),r)=pDoir) [1 - ap(r)] , 



(11) 



where a > is the first order correction coefficient. In 
(|ll|l we assume that the HI between A and B, expressed 
in Do{r), are independent of the HI between B particles 
and the resulting mobility can be factorized. For the 
case of hard spheres with diameter a it is well established 
that aa~^ ~ 3.43»i^ Substituting Ullfl and the derivative 
of with respect to p into Eq. ©, integrating and 
strictly linearizing in po we find 



C 



Dairy 



p{r) 



l3dVAB{r) 
dr 



2BXC 



Doiry 



d_ 

dr 



pM,(12) 



where we introduced the effective second virial coefficient 



B* 



B2 - a/2. 



(13) 



and which has to be solved with boundary conditions 
Eqs. (|7|) and (jSJ. The equilibrium pair interaction, de- 
scribed by B2, and the dynamic interaction, absorbed 
in a, compete likewise now in altering the rate of the 
reaction, and can hence be summarized in the single pa- 
rameter B2 given in (|13l) . This interesting feature will be 
discussed in more detail in the next section, where the 
result for k is presented. Eq. (|12|l is a first order lin- 
ear differential equation and can be solved analytically. 
The parameter C is an integration constant, determined 
by the boundary conditions, and is related to the rate 
constant with C = pok/An. 



B. Result for the steady-state rate constant 

From the linearized solution of 1)12(1 we arrive at the 
result for the diffusion-controlled steady-state rate con- 
stant valid for small densities (or weak interactions) 



1 + B^po 



/i(^) 



(14) 



where fc^i is the classical result of Debye^ for ideal or 
infinitely diluted substrates in a nonzero potential Vab (r) 



k 



D 



dr 



and 



/(r) = 2 



47rL'o(r)r2 

Io{r)hir) ~ hir) 



/i(oo) 



{/i(^)}-i (15) 



(16) 



with the integrals 



1 



h{r) 



dr' 



4TTDa{r')r'^ 



and 



hir) 



dr' 



^TTDoir'y^ 



(17) 



(18) 



(19) 



which are functions depending on the particular inter- 
action VAsir) and the position-dependent diffusion con- 
stant Do{r). When interactions between A and B can be 
neglected, namely Do{r) = Do, and Vab(?') = 0, it fol- 
lows /(cxd)//i(oo) = 1, and the result for the steady-state 
rate constant reduces to 



k = A7rDo(TABil + B^pq), 



(20) 



which for noninteracting particles or at infinite dilution 
of species B gives the classical Smoluchowski result^ 



fco = AttDoctab- 



(21) 



The constant /(oo)//i(oo) in l|14|l depends on the par- 
ticular interaction Vab(?') and Df){r) between A and B 
particles, but is always close to 1 for Lennard- Jones or 
Yukawa like interactions typically found in solutions, and 
a Do{r) for instance given by the Oseen tensor.^ The ma- 
jor contribution to the density correction stems from the 
effective second virial coefficient (|13f) , which is a corrected 
second virial coefficient for the dynamic situation and fea- 
tures the following interesting trends for the steady-state 
rate constant. Generally a > 0, imposing that HI will de- 
crease the substrate mobility. However, if the substrate- 
substrate interaction is strongly repulsive, B2 > a/2, so 
that i?2 is positive, an enhanced reaction rate is predicted 
compared to the noninteracting case. On the contrary, 
a mainly attractive pair interaction will in general lead 
to a negative B2 and decrease the reaction rate. For 
long-ranged interactions |i?2| may be an order of magni- 
tude larger than a/2 and HI can be neglected, such that 
B2 — i?2- When the repulsive interaction leads to a B2 
comparable to a/2, both contributions can cancel each 
other and the reaction rate just changes little or not at all 
with increasing (small) density or (weak) interaction. For 
instance, for the repulsive hard sphere case the effective 
second virial coefficient B^cr'^ ~ 2.09 - 3.43/2 ~ 0.38 is 
small due to a slowing of the diffusion of the substrates 
which is comparable to the effects of the enhanced os- 
motic pressure in the system, leading just to a minor 
increase of the rate constant with density. Our theory 
is in agreement with previous workiSii^ where the ex- 
cluded volume of the substrates is predicted to lead to an 
enhanced reaction rate (excluding HI). It is worth men- 
tioning that B2 is half of the coefficient of the linear first 
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order correction in po of the collective diffusion coefficient 
of the substrates in bulk solution, latter evident from the 
generalized Stokes-Einstein relation ||3Jl. The range of 
validity of H14(l will be examined in section III for hard 
sphere and repulsive and attractive Yukawa interactions 
excluding HI. 



C. Result for the density profile 

The general solution for the diffusion-controlled 
steady-state density profile in linear order in po reads 



P{r) = Po 



kMr) 
exp[/3VAB(? 



1 + 



I{oo) I{r) 



/i(oo) /i(r) 



(22) 



In order to better identify how the substrate interac- 
tions affect the density profiles we rewrite the solution for 
the case of no interactions between substrate and sink, 
Do{r)=Do and Fab (r) =0: 



where 



p{r) = Pid (1 + Bjpo^^) , 
iir) = PO (1 



Pid I 



CAB 



(23) 



(24) 



is the classical result for ideal or infinitely diluted sub- 
strate particles when Vab = and DQ{r) — Dq"^ Eq. 
H23|l predicts that a positive B2 will enhance the sub- 
strate density close to the absorbent due to a higher 
(positive) bulk pressure of the B particles compared to 
the noninteracting case. This is opposite to the case of 
negative B2 (mainly attractive interactions) where a de- 
pletion of B particles takes place close to A due to a 
negative bulk pressure, and/or a larger immobility of the 
approaching substrates. Comparing these trends to the 
rate constant l(T^ we conclude that an enhanced or lower 
density close to the sink increases or decreases the reac- 
tion rate, respectively, at least for small densities. These 
predictions will be examined in the next section with BD 
simulations. 



III. BROWNIAN DYNAMICS SIMULATION 

A. Simulation details 

For a verification of the theory in the previous sec- 
tion we perform standard Brownian Dynamics simu- 
lations using the integration technique of Ermak and 
McCammon;^ and in addition accounting for hard 
sphere overlap^ In order to conduct a clear and trans- 
parent comparison to the theoretical result we neglect 
HI in the BD simulations. An accurate treatment of HI 
is only given by more sophisticated and computation- 
ally expensive techniques, such as Lattice-Boltzmann^*^ 
or others iSiiSSiS^ Thus, in the following Do{r) — Dq given 



by and i?2 = S2. In the BD simulations, the steady- 
state rate constant is determined by calculating the sur- 
vival probability S{t) of the sink particle A in the steady- 
state, which is given by 



S{t) — exp(— fcpo^) 



(25) 



in the presence of reactive B particles with bulk density 
Po- In our simulation model, the sink is fixed in the center 
of the periodically repeated, cubic simulation box with 
length L. Substrate (B) particles react with the sink (A) 
as soon they touch it, which happens at a center-to-center 
distance tab — cab- The trajectory of a B particle is ter- 
minated after reaction. The survival probability can be 
easily measured by forming an histogram of the times be- 
tween two reactions and normalize it with respect to the 
number of total reactions in one simulation. To account 
for the steady-state case and simulate a fixed density po 
far away from the sink, the annihilated B particle is rein- 
serted randomly at a position close (rAB > 0.45L) to 
the box edges. Due to the long range of the steady-state 
density-profile H23I) , finite size effects can be large for too 
small box sizes for all densities and have been analyzed 
by performing finite size scaling simulations. Eventu- 
ally the box lengths used in our simulations range from 
L = 50(7 up to L = 200(7 depending on density and in- 
teraction of the B particles. This involved simulations of 
Nb = 1000 - 30000 substrate particles using 5 x lO"^ to 
2 X 10^ timesteps. The integration timestep was chosen 
to be 0.003ts, where tb — cr'^/Do is a typical Brownian 
timescale in the simulation. Verlet-neighbor lists were 
used to optimize computational time.^§ Densities up to 
B2P0 ^ 1 could be simulated; for larger densities the 
system size and statistical errors become too large for 
reasonable output. 



B. Systems 

In order to perform a systematic comparison to the 
theoretical results we consider four different systems, I, 
II, HI, and IV which differ in the interaction between 
A and B, and B and B particles. In all systems the 
interactions are in general given by 

PV,, (r) - pV^s [Uj ) + U,j ^ exp[-K(ry - )] (26) 



where /3Vhs is the hard sphere interaction 



for < CTy ; 



else, 



(27) 



and the second term in H2t)|) is a Yukawa interaction 
with energy scale Uij, inverse screening length k, and 
i,j = A^B. K will be fixed to kct = 1 in all our four 
systems. We also fix the length scale to ctab = l-Scr, 
which is given when the sink particle is twice as large as 
the substrate particle ctaa = 2(7. We note already that 
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System 


f/AB 


f/BB 


o"ab/o" 


B2cr-3 


7(oo)//i(oo) 


I 


0.0 


0.0 


1.5 


2.09 


1.00 


II 


0.0 


1.0 


1.5 


13.32 


1.00 


III 


-1.0 


1.0 


1.5 


13.32 


1.11 


IV 


0.0 


-0.5 


1.5 


-4.62 


1.00 



TABLE I: 



we have simulated selected densities for all four systems 
for two other sink sizes, namely cab — c and ctab — 2cr, 
and did not find any deviation from the results later in 
the work. For high asymmetries (Taa ^ which is the 
more realistic case, the BD simulations unfortunately be- 
come computationally expensive due too a large number 
of substrate particles which have to be simulated. Note 
again that the range of the steady-state density profile 
increases essentially linearly with ctab, see Eq. (|^ : how- 
ever, we stress that this does not affect our theoretical 
result H14() which should be valid for all size ratios, and 
we have chosen smaller sink sizes only for computational 
convenience. 

In system I we simply consider Vij = 0, in which case 
all interactions are hard sphere like. The second virial 
coefficient for HS is B2<7~^ = 27r/3 ~ 2.09. In system 

II Yukawa-like repulsion C/bb = 1 is added between the 
substrate particles. The second virial coefficient increases 
to a total of S2CT-3 ~ 27r/3 -f 11.23 ~ 13.32. In system 

III a Yukawa attraction is added between sink (A) and 
substrate (B) particles, Dab = —1, while keeping the 
B-B interaction as in system II. Finally, in system IV 
we choose J/ab = 0, no interaction between A and B, 
but here we focus on attraction between the B particles, 
C7bb = —0.5, resulting in a total negative second virial 
coefficient B2(J~^ — —4.62. The system parameters and 
corresponding values of the second virial coefficient and 
the constant I{oo)/ Ii{oo) are summarized in Tab. 1. 

Typical examples of the calculated survival probabili- 
ties (|25ll are shown in Fig. 1 for systems I, II, and IV on 
a logarithmic ordinate and show linear behavior for all 
times within the uncertainty of the statistics. The den- 
sity for system I (hard spheres) is chosen to be very small 
{B2P0 = 0.001) and the survival probability matches with 
the classical theoretical result S{t) = exp(— fcopo^), also 
shown in Fig.l as dashed line. As expected from the 
theory, a positive or negative B2 (systems II and IV) in- 
creases or decreases the rate constant, respectively, as is 
indicated by an increased or decreased absolute slope of 
the S{t) curves. Linear fits through these curves deter- 
mine our 'experimental' values of the rate constant for 
all systems, while the regression coefficient of the fitting 
procedure provides the error bars. 




FIG. 1: Typical examples for the survival probability S(t) on a 
logarithmic scale versus time scaled by the ideal substrate rate 
constant ko 1211 1 and substrate bulk density po- The straight 
dashed line is the ideal result S{t) — exp(— fco/Oot), while the 
noisy data are BD simulation results for chosen densities in 
systems I, II, and IV. 



C. Results 

Examples of the density profiles for system I, II, and IV 
are shown in Fig. 2 compared to the theoretical prediction 
(|23|l and are in very good agreement. The density for sys- 
tem I is in the very dilute region {B2P0 — 0.001) and the 
profile matches the 1 — ctab / ''ab behavior H24|l for ideal 
substrate particles. Further on, the predicted density 
increase (decrease) around the sink by repulsive (attrac- 
tive) substrate interactions is verified by the BD simula- 
tions. All observed density profiles exhibit a very small 
but nonvanishing density value at contact /9(ctab) — 0.05, 
showing that the assumption p{a ab) — for absorbing 
boundary conditions in the theory is justified, but not 
perfect. This contact value increases (decreases) slightly 
with positive (negative) 

Results for the rate constant are shown in Fig. 3 where 
we plot the steady-state rate constant scaled by the ideal 
rate constant fcg given by H21|l versus the substrate den- 
sity scaled by the second virial coefficient B2, which can 
be different for each system, see Tab. 1. Note again 
that we neglect HI in our simulations, such that a = 
and Do{r) = Dq, and B2 = B2 in the theory. We find 
excellent agreement to the theory (lines) within the sta- 
tistical uncertainties of the simulation for the systems 
with repulsive substrate interaction (I-III) up to densities 
B2P0 — 0.4. For system I, the hard sphere case, the non- 
linear regime takes over for larger densities and the the- 
oretical prediction underestimates the rate constant. For 
systems II and III we still find good agreement (theory 
within the error bars) till up to B2P0 — 0.7. This agree- 
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FIG. 2: Steady-state density profiles p(r) of the substrates 
around the sink particle from BD simulations (symbols) and 
theory (lines) according to Eq. 12311 for three different pa- 
rameter sets, r is the distance from the center of the sink. 
The circles are close to the infinite dilution limit l|24^ . while 
squares and diamonds are for positive and negative B2 for 
nonzero bulk densities po, respectively. 



tem with attractive substrate interactions (IV) the agree- 
ment is good only for smaller densities {B2P0 < 0.20). 
EquiUbrium statistics has shown that for attractive inter- 
actions, a LDA approach usually fails to give an accurate 
description of the system behavior even for weak interac- 
tions, which could also be the case here. Secondly, finite 
size effects in the simulation are larger due to a longer 
ranged density profile (|23(l and cannot be excluded to 
explain the discrepancies for larger densities. The theo- 
retical curves for I and II are on top of each other and 
have the same absolute slope as the data for IV because 
the density is scaled by the corresponding For zero 
density the curves for I, II, and IV intersect at fc/fco = 1 
which is the classical ideal gas limit for no interactions 
between B and B, and A and B. Curve III intersects at a 
larger value k/ko ~ 1.26 due to the additional A-B inter- 
action in this system. This is also the main reason for the 
increased absolute slope compared to I, II, and IV, but 
additionally the constant /(oo)//i(oo) ~ 1.1 is enhanced 
by the attractive A-B interaction in this system. 



IV. CONCLUDING REMARKS 




0.25 0.5 0.75 



FIG. 3: Dimensionless reaction rate k/ko from BD simulations 
(symbols) and theory (lines. Eg. 11411 1 versus density of the 
substrate particles for the four different systems I (circles) , II 
(squares). III (diamonds), and IV (triangles). Note that the 
density is scaled by the corresponding second virial coefficient, 
see Tab. I, that is why the theoretical curves for I and II are 
on top of each other and have the same absolute slope as the 
curve for IV. 



ment for a larger density range as compared to HS could 
be anticipated since for soft repulsive interactions it is 
known that higher order virial coefficients (-B3, and so on) 
are less important even for higher densities;^ However, 
for larger densities the shortcomings of the linearized 
Smoluchowski equation also take their effects. In the sys- 



In conclusion we have derived an analytic expression 
for the density profile and rate constant for weakly in- 
teracting substrate particles for the steady-state case of 
diffusion controlled reactions. For this purpose we used 
the Smoluchowski equation, which was generalized within 
a LDA to account for the osmotic pressure and HI of the 
interacting substrates. A comparison to BD simulations 
excluding hydrodynamic interactions (accounted for in 
the theory by a = 0) showed excellent agreement for 
densities up to B2P0 — 0.4 for repulsive interactions, and 
up to B2P0 — 0.20 for Yukawa-like attractions. Our BD 
simulations do not include HI but support our theoretical 
concept and prove its validity when HI can be neglected, 
i.e. for long-ranged interactions. We are confident that 
our theoretical treatment is valid also in the general case 
including HI; when the density-profiles of the substrates 
are slowly varying in space, the correlations in the sys- 
tem are usually well approximated by those in the bulk, 
even in the dynamic case.^^ However, a verification of our 
theory including HI is highly desired and abandoned to 
future work, where the HI need to be treated by accurate 
means. 

In principle, our study provides an alternative way 
to estimate the second virial coefficient of interact- 
ing macromolecules experimentally by measuring their 
steady-state rate constant in diffusion-controlled reac- 
tions at low densities. In such experiments, the effec- 
tive sink-substrate interaction has to be known (e.g. by 
measuring the infinite dilution limit of the rate constant) 
and the hydrodynamic quantities Do{r) and a must be 
approximated, e.g. by using the Oseen tensorj^ and ap- 
proximating the substrate particles by hard spheres with 
an effective hydrodynamic radius, respectively»i& 

The generalized Smoluchowski equation used in our 



7 



work can also be understood as a LDA in the recently 
proposed framework of dynamic density functional the- 
ory (DDFT),-^ where equilibrium correlations are used 
to approximate the dynamical correlation in a Brow- 
nian system. It was shown by BD simulations that 
DDFT including more sophisticated approximations than 
LDA works well in the case of dense one28 and three- 
dimensional hard sphereS)2S and three-dimensional par- 
ticles with very soft interactions )22i2i*^ and seems to 
provide a powerful tool to extend our work to more 
strongly correlated systems, i.e nucleating or aggregat- 
ing colloids, polymerization, or binding in crowded pro- 
tein solutions 1^ 

Finally, we hope that our approach will shed some 
light on interpretations of experimental rate con- 
stant measurementSi^ and will be useful in extend- 
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ing existing works on the time-dependence of the rate 
constant crowding effects in reactions ji2*ii finite 
sink concentrationj^ii^ii^ or anisotropic reactivitji^i^ to 
the case of interacting substrate particles. 
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